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Abstract. We characterize the multipartite entanglement of a system of n qubits in 
terms of the distribution function of the bipartite purity over all balanced bipartitions. 
We search for those (maximally multipartite entangled) states whose purity is minimum 
for all bipartitions and recast this optimization problem into a problem of statistical 
mechanics. 
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Entanglement is a consequence of the quantal superposition principle and therefore 
of the linearity of Hilbert spaces, and embodies the impossibility of factorizing a given 
state of a total quantum system in terms of the states of its constituents. The notion of 
entanglement plays a central role in the development of novel resources and techniques in 
computing, communication and information processing. There is a widespread literature 
about the characterization, quantification and properties of bipartite entanglement, i.e. 
the entanglement of two subsystems [U[2]. The landscape of multipartite entanglement 
is less understood but very widely investigated [31 HI El E] . The main difficulty is that 
there is no unique way of quantifying it: more so, different definitions often do not agree 
with each other, because they adopt different strategies and capture different features 
of this inherently quantum phenomenon. 

It would be desirable to understand how to properly characterize multipartite 
entanglement, e.g. by identifying a few key properties that can account for its overall 
features. The quantification and evaluation of global entanglement measures bears 
serious computational difficulties, because states endowed with large entanglement 
typically involve exponentially many coefficients. Multipartite entanglement has 
therefore the typical features of a complex phenomenon, and for this reason we shall 
follow a precept derived from the study of complex systems, that relies on the idea 
that complicated phenomena cannot be summarized by a single (or a few) number(s) 
[7j. We notice that in the context of quantum entanglement this idea was alluded to in 
[8]. In this spirit, entanglement can be analyzed in terms of a function: the probability 
density of the purity of a subsystem over all bipartitions of the total system [HI [TO] . 
A state has a large multipartite entanglement if it has a large and well distributed 
bipartite entanglement. In other words, states with large multipartite entanglement are 
characterized by a narrow distribution of bipartite entanglement, centered around a 
large average. 

In this letter we shall explore these ideas even further and endeavor to give a 
statistical mechanical characterization of multipartite entanglement, by recasting the 
problem in terms of a cost function that will represent the (inverse) global entanglement 
of a quantum system. We will establish a direct connection between the study of 
multipartite entanglement and the tools of classical statistical mechanics. 

We shall consider an ensemble of n spin- 1/2 particles (qubits), whose Hilbert space 
% = (C 2 )® n has dimension N = 2 n . We divide this set in two parts, A and A, made up 
of ua and qubits respectively {ua + nj — n). The total Hilbert space is split in the 
tensor product H = Ha ® Hli with dimensions Na = 2™ A and Na = 2"^ respectively 
(NaNa = N). We also assume with no loss of generality ua ^ ua- We shall consider 
only pure states and shall not discuss additional phenomena such as bound entanglement 
[m[T2]. Let 

N-1 N-1 

hA> = 5>|j>, Zj eC, 5>/ = i, (l) 

j=0 j=0 

where the computational base {\j)} is expressed in terms of the eigenstates of the third 
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Pauli matrices acting on each site. These yield binary sequences belonging to {0, l} n , 
whose decimal representation are the indices j. For each bipartition (A, A) we define 
a bijection j -H- (j'a, Ja); where the integers ja G [0,N a — 1] and ja G [0, Na — 1] are 
the decimal representation of the subsequences of portion A and A, respectively. This 
reflects the factorization \j) = \j A ) <8> \Ja)- As bipartite entanglement measure [at fixed 
bipartition (A, A)] we take the purity of subsystem A, which is conveniently expressed 
in terms of the reduced density operators, pa = tr^ p (p = \ip)(ip\ being the total density 
operator) 

7l A = tr A p 2 A = tr A PA = Z 3A,j A ^A,3A Z lA,lA^A,lA- ( 2 ) 

3A,Ia ,3A' l A 

Purity ranges in [1/N A , 1] and is 1 if and only if pa and pa are projectors, that is, if 
the state is (bi) separable, \ip) = \i/ja) <8> \4>a)- Otherwise \ip) is entangled, tia saturates 
its minimum 1/Na if and only if the reduced density operator of the smaller partition 
is proportional to the identity matrix pa = 1/N A - The latter is the case of states 

endowed with maximum bipartite entanglement, at fixed bipartition (A, A). The 
distribution function of purity over all bipartitions, p(tta), yields a probability-density- 
function characterization of multipartite entanglement: its average measures the mean 
entanglement of the state when the bipartitions are varied, while its variance quantifies 
how uniformly bipartite entanglement is distributed among all possible bipartitions. 

The coefficients Zj in ([1]) belong to the set 

C = |(^ , zi, ... , ztf-i) e C N \ N 2 = 1 J ■ (3) 

Note that C corresponds to the set of normalized vectors 5 = {^G'H|||'?/ ; || = 1} and 
is left invariant under the natural action of the unitary group tl(T-L). A typical state is 
obtained by a uniform sampling of S. Typical states can be (efficiently) generated by a 
chaotic dynamics [TBI fUj . 

In the limit of large N, the vr^'s have a bell shaped distribution over the bipartitions 
with mean and variance 

Na + N a 

A* = (*a) = N+1 , (4) 

2 /t ,2, 2{N\-l){N\-l) 

a =((*A-H))= {N+mN + 2){N + 3 y (5) 

respectively jSHinilI51IISlinilIHllISll2D]- The brackets (■ • •) denote the average with 
respect to the unitarily invariant measure over pure states 

dpc(z) = Y[ dz k dz k 5 j 1 - ^2 \ z k? ) , (6) 

k \ k J 



induced by the Haar measure over U(H) through the mapping \ip) = Zj\j) = U\i(jo), 
for a given state |^o) [18] . 

It is interesting to notice that, for balanced bipartitions ha = /2 J , the average 
and standard deviation of the distribution simplify to pt, ~ 2/Na,ct ~ y/2/N\, the 
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Figure 1. Probability density function p(n A ) of a typical state over balanced 
bipartitions. For large n, the average is 2/Na and the standard deviation \/2/N\ 



{Na — \N). Our objective is to characterize those "maximally multipartite entangled 
states" whose average purity is minimum. 



approximation becoming quickly very accurate as the number of qubits increases. The 
behavior of the probability density function p(tt a ) is sketched in Fig. [TJ It is worth 
emphasizing that the typical states are very entangled, but are far from the minimum 
(on average) by a factor 2; moreover, their distribution function is very narrow (a 
signature of a good multipartition of entanglement), but has a nonvanishing standard 
deviation. It is then natural to ask: Is it possible to find states that perform better? The 
answer to this question is positive [21]: there exist "maximally multipartite entangled 
states" (MMES) such that their distribution function (over balanced bipartitions) is a 
Dirac delta function. For n = 2, 3, 5, 6 the delta is centered on the smallest possible 
value: p(tt a ) = 5(tt a — F° r n = 2 MMES are Bell states up to local unitary 

transformations, while for n = 3 MMES are equivalent to the GHZ states [22]. For 
n = 4 one numerically obtains p(n A ) = 5(ii A — C 4 ), with C 4 = 1/3 > 1/4 = 1/Na 
(but still C A < 1/2 = 2/N A ) [2D EU E3 [25]. For n = 7 one finds states such that 
p(7r A ) = 5(tt a - C 7 ), with C 7 ~ 0.136 > 1/8 = l/N A . "Perfect" MMES (that saturate 
the minimum purity 1/N A ) do not exist for n > 8 [26]. This will be presently viewed 
as a symptom of frustration and is a new and intriguing feature of the multipartite 
scenario, as opposed to the bipartite case. We shall now endeavor to outline a strategy 
that enables one to characterize MMES. 

The search for MMES can be recast in terms of an optimization problem. Look at 
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the average purity as a cost function ("potential of multipartite entanglement" |2T|) 



where E denotes the expectation value and the sum is over balanced bipartitions 
ha = |_ n /2_|- Since we are focusing on balanced bipartitions, and any bipartition can 
be brought into any other bipartition by applying a permutation of the qubits, the sum 
over balanced bipartitions in (ED) is equivalent to a sum over the permutations of the 
qubits. H in Eq. (J7J) is related to the linear entropy considered in (261 127] . 

The above minimization problem has proved to be of formidable difficulty, mainly 
for the presence of frustration: the requirement that purity saturates its minimum 1/Na 
can be satisfied for some but not for all balanced bipartitions. In other words, not all 
terms in the summation (JTj) can be 1/Na- In order to study this problem, we will recast 
it in classical statistical mechanical terms: instead of searching for the minimum of H, 
one can look at the free energy of a suitable classical system at a fictitious temperature 
and recover the original problem in the zero temperature limit. 

We follow the standard ensemble approach of statistical mechanics and introduce 
an ensemble {mj} of M vectors (states), where mj is the number of vectors with purity 
H = €j. One seeks the distribution that maximizes the number of states Q = Ml / Y[j m 3 -! 
under the constraints that Ylj m j = M and Y2j m j e j = ME. For M — > oo, the above 
optimization problem yields the canonical ensemble and its partition function 



where the Lagrange multiplier j3, that plays the role of an inverse temperature, fixes 
the average value of purity E. The constraint C in the integration domain is due to the 
normalization ([3]) and the measure is defined in ([6]). Note that the partition function 
([8]) can be given a base-independent expression 



Z{/3) = / al^c{z)e- mz) = c N / dfi n (U)exp (- (3E[tr A (tr A U\^)(MU ] f]) , (9) 



where //h denotes the Haar measure over lf(/H), |^ ) is an Y given vector and the 
(unimportant) constant cn = /ic(C 7V )//iH(W('H)) is proportional to the ratio between 
the area of the (N — l)-dimensional sphere ([3]) and the volume of the unitary group. 

The average over bipartitions in (0) entails a complicated dependence on the indices 
j in Eqs. (CQ)- (TSJ). This hinders one from finding a closed expression for the partition 
function and makes the application of a standard saddle-point method or the calculation 
of the quenched average rather involved. Nonetheless, by applying the standard tools 
of statistical mechanics, one can discuss several aspects and analyze interesting limits. 
For (3 — > Eq. ([8]) yields the distribution of the typical states (Q. Notice that this 
is valid both for (3 — > O 1 * 1 . For /3 — > +oo, only those configurations that minimize the 
Hamiltonian survive, namely the MMES. Remarkably, there is a physically appealing 
interpretation even for negative temperatures: for — > — oo, those configurations are 




(7) 




(8) 
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Table 1. High, low and negative temperature limits. 



13 -> +oo 


# = 


E'o (min) 


MMES 


p -> 






typical states 


fi ->■ — oo 


# = 


1 (max) 


separable states 



selected that maximize the Hamiltonian, that is separable states. These preliminary 
findings are summarized in Table [TJ 

By manipulating the partition function we can write the energy distribution 
function at arbitrary fi in terms of its high temperature limit: 

PJE) = — l-rr f diJc(z)5 (H - E) e~ m = 6 (io) 

; Z(J3)J ' 1 ; jldEe-P E P (E) 1 ; 

where G [i?o> 1]> -^o being the minimum of the potential of multipartite entanglement 
H in (J7J). We know that 1/JVa < S (JV A ) < fi < 2/N A , and thus lim^^oo E Q (N A ) = . 
By multiplying and dividing the last equation by \(i\e^ and fie 13 ^, respectively, we find 

P- 00 (E) = 5(E-1), P +oo (E) = 5(E-E ). (11) 

These limits are the counterparts of those discussed for the partition function and are 
reflected in the asymptotic behaviour of the average energy as function of fi 

= Y^)j d^c{z)He- pH = jf alEEP p {E) = ~\*Z{p). (12) 

Indeed, 

(#)_«, = 1, (H) +00 = E . (13) 

Moreover -§p(H)p = ~{H 2 }p + (H} 2 p = -AHj ^ 0, which is nonpositive. Thus the 
average energy is a non-increasing function of fi and has at least one inflexion point 
as function of fi. From a qualitative point of view, one expects the behavior sketched 
in Fig. [2) for j3 — > 0, the distribution is Gaussian (typical states); when fi — » +oo 
the distribution tends to become more concentrated around E . This picture will be 
substantiated in the following. 

In general, the high-temperature expansion of the average energy reads 

°° f o\m—l 

m=l V >' 
am— I 

K M {H) = { _ r -i^__ {Hhj (15) 

where ^^(H) are the cumulants of H. So far, the analysis is valid at fixed N = 2™. 
Let us now briefly discuss the limit iV — > oo. The first cumulant is nothing but the 
average purity [i defined in (j3J). The second cumulant a 2 = AHq = Kq(H) can also be 
computed exactly. A detailed calculation, to be presented elsewhere, yields 

a 2 = 4 2) (#) ~ 3V2N~ 4+1 °^ 3 ~ 0{N~ 2A2 ). (16) 
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Figure 2. Qualitative sketch of Eq. (|T0|) . at fixed Na, in arbitrary units. The energy 
density function is distributed around /i with standard deviation a at j3 = (inset) 
and moves toward E when (3 — > +oo. From right to left, j3 changes in constant steps. 
The probability density rigidly shifts with /3, for < N 7 / 2 - l °^ 3 . See Eq. JIS) and 
following discussion. Note that E Q = 0(N~ 1 / 2 ). 



The presence of an irrational exponent is a rather rare feature and an interesting 
phenomenon in itself. We notice that if the vr^'s in the sum (J7|) were independent 
Gaussian random variables, one would get a 2 ~ a 2 /N = 0(N~ 3 ), where a is defined in 
Eq. (jSJ) (the number of balanced bipartitions being of order N.) This shows that different 
bipartitions "interfere," yielding a highly nontrivial average. In fact, it can be shown 
that the correlations decrease exponentially with the distance between bipartitions. 

Higher order cumulants at (3 = can be shown to tend to zero faster than the others. 
Since the Gaussian is the only probability distribution with a finite number of non- 
vanishing cumulants (the first and the second one), one can therefore assume a Gaussian 
approximation for the energy distribution at infinite temperature, by neglecting all the 
cumulants but the first two: 



where \i and a are defined in (@J and (TT6|) . respectively. Finally, by applying ( 1TU1) we 
can compute the energy distribution at arbitrary temperature 



Thus from the definition (fT5l) of cumulants, {H)p as a function of (3 would have a 
vanishing curvature, in contradiction with the qualitative sketch mentioned after Eq. 




(17) 




(18) 
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(fT2j) . Eq. (USD can be written as 

P l3 (E) = P (E + f3a 2 ), (19) 

an expression valid also for nongaussian distributions that are well concentrated around 
their mean /i. As a matter of fact, Eqs. (fl8l) and (fT9l) are valid only when [3 is not too 
large. In fact, they hold true as far as /i — (3a 2 — a > 0, i.e. (3 < /i/(x 2 = 0(iV 7 / 2 ~ log23 ) ~ 
0(iV L92 ). Up to this value the probability density rigidly shifts with /3, as is apparent in 
Fig. [2j For larger values of (3 the lower tail of the distribution starts "feeling" the wall 
at Eq. The large-/? asymptotic form of Pp(E) depends on the behavior of Pq(E — > E ): 
one easily obtains 

ai+i 

Pp(E) ~ —^-(E - EoYe-W-**, (20) 

where ^ is the order of the first nonvanishing derivative of Pq(E) at Eq. (Figure [2] 
displays the case t = 0.) Notice that the only relic of Pq{E) in (j20|) is £ and P^ +OQ (E) 
yields the second equation in (TTTT) . The analysis for /3 — > — oo is analogous and yields 
the first equation in (fTTj) . 

In conclusion, we have built a statistical theory around an optimization problem, 
whose solutions are the maximally multipartite entangled states, that appear as minimal 
energy configurations. The approach we propose borrows methods from classical 
statistical mechanics in order to investigate the multipartite entanglement scenario. 
The introduction of a temperature is a familiar expedient in statistical mechanics, for 
instance in the study of optimization problems with simulated annealing and tempering. 
One obtains here an interesting picture for all real values of (3 that fixes, with an 
uncertainty that becomes smaller for larger systems, the value of the purity of the 
subsystem under consideration, thus identifying an "isoentangled" submanifold of states. 
A strategy similar to the one adopted here was used in [28] for the simpler case 
of bipartite entanglement at a fixed bipartition, where the purity exhibits a phase 
transition. We have seen that the multipartite version of the problem is much more 
complicated, but it would be of great interest to understand whether the phase transition 
that occurs in the bipartite situation, when there is no average over the bipartitions, 
has a counterpart in the multipartite scenario. This would help elucidate the relation 
between multipartite entanglement, complexity and the presence of frustration. 
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